Superconductivity in correlated disordered two-dimensional electron gas 
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We calculate the dynamic effective electron-electron interaction potential for a low density dis- 
ordered two-dimensional electron gas. The disordered response function is used to calculate the 
effective potential where the scattering rate is taken from typical mobilities from recent experi- 
ments. We investigate the development of an effective attractive pair potential for both disordered 
and disorder free systems with correlations determined from existing numerical simulation data. 
The effect of disorder and correlations on the superconducting critical temperature T c is discussed. 

PACS numbers: 73.20.Dx,71.30.+h,71.55.-i 



Recent experiments have revealed a well denned metal- 
insulator transition (MIT) in 2D electron systems. The 
MIT has been observed in number of different systems, Si 
MOSFET §, GaAs/AlGaAs {§ and p-SiGe §) for both 
electrons and holes. A common feature near the transi- 
tion is that the average electron-electron interaction en- 
ergy is an order of magnitude larger than the Fermi en- 
ergy indicating that the transition is driven by the strong 
Coulomb interactions. 

The transition from the metallic phase has been found 
to occur at widely different levels of disorder and mo- 
bility. In Si MOSFETs the disorder at the transition is 
strong [Q while in high mobility SiGe samples H the 
disorder is very small. 

Another feature of the MIT is that the resistivity scales 
with temperature § near the transition according to 
p(T,n s ) = p(T/T ). The single parameter T a is a func- 
tion of 8 n = (n s — n c ) /n c on both the metallic and insu- 
lating sides of the transition. n s is the electron density 
and n c the critical density at the transition. As a re- 
sult of this scaling all the data can be collapsed into two 
branches, an upper one for the metallic phase and a lower 
one for the insulator. The scaling behavior of the resistiv- 
ity is also observed with external electric field |J. Such 
scaling is characteristic of a true phase transition. It is 
similar to the superconductor-insulator phase transition 
in thin disordered films and the quantum Hall liquid- 
insulator phase transition. Independent measurements of 
the I-V characteristics at low temperature and the mag- 
netic susceptibility are required before the existence of a 
superconducting phase can be confirmed. 

Thakur and Neilson Q have proposed that the insulat- 
ing phase of the low density system with weak disorder 
would be a frozen solid phase with liquid-like short-range 
order. They showed such a phase could exist for electron 
densities corresponding to r s >7 at impurity levels which 
were in good agreement with the observations of the MIT. 

The metallic phase near the transition is not a con- 
ventional metallic state for which the Coulomb interac- 
tion energy is quite small compared to Fermi energy. If 
that had been the case properties could have been deter- 
mined perturbatively in powers of the interaction poten- 
tial. However, near the MIT the electrons are strongly 



correlated on both the metallic and insulating sides. The 
scaling of the conductivity on both sides with a common 
parameter T Q indicates that the mechanisms for electri- 
cal transport in the insulating and metallic phases are 
related. 

A number of proposals have been made about the na- 
ture of the metallic phase. One suggestion is supercon- 
ducting pairing of the electrons with the pairing me- 
diated by the dynamic correlation hole surrounding each 
electron. At very low densities the correlation hole and 
the electron move rigidly as a single unit, and in the ex- 
treme case of the Wigner crystal the correlation hole and 
the electron form the Wigner-Seitz cell. In the conduct- 
ing phase at comparatively higher densities than this the 
correlation hole is not rigidly bound to its electron and 
on a time scale of the inverse of the plasmon frequency 
there is the possibility of relative movement leading to 
partial dissociation. Over this time scale the partially 
vacated correlation hole can attract another electron. 
This attraction can lead to superconducting pairing of 
the electrons. The existence of a strong correlation hole 
and weakly damped plasmons are both necessary for this 
mechanism. When the correlation hole is weak or disor- 
der strongly damps the plasmon there is no attraction. 

In this paper we investigate the combined effect of dis- 
order and correlations on the superconductivity. We find 
that the superconductivity persists at levels of disorder 
where the plasmon is damped out. We also find that the 
superconducting transition temperature is more sensitive 
to disorder when the electron correlations are weaker, 
that is at higher electron densities. 

To treat the disorder we introduce the generalized ran- 
dom phase response function 



1 + V(q) [l-G(g)] X « (<?,«/) 



(1) 



The local field factor G(q) takes into account the corre- 
lations. These make the effective electron-electron inter- 
action weaker than the bare Coulomb potential V(q) = 
2ne 2 /qe. The effect of disorder is contained in the suscep- 
tibility (<7, oS). This is the particle conserving response 
function of non-interacting electrons scattering from the 
disorder O], 
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7 is the scattering rate off the disorder. X (l) is the 
static susceptibility for non-interacting electrons with- 
out disorder. x fai^ + ry) is the dynamical suscep- 
tibility for non-interacting electrons scattering off the 
disorder. In the diffusive regime lim Wj9 ^o X (l, w ) — 
{2m*/irk F h 2 ) {Dq 2 /(Vq 2 + iw)}, where V = v 2 F /~f is the 
diffusion constant. This vanishes as 7 goes to infinity, 
representing the transition to the insulating phase ]T2| . 

The scattering rate 7 could be calculated using the 
memory function formalism for the density relaxation 
function and applying mode coupling theory [ O] . How- 
ever, here we take 7 as a parameter, the value of which 
is representative of the experimental mobilities observed 
in recent studies. 7 is related to the electron mobility /1 
by the Drude expression /1 = e/m*j. 

The local field factor G(q), taking into account electron 
correlations, is calculated from numerical simulation data 
of the static structure factor S(q) Qj. Introducing the 
fluctuation dissipation theorem, G{q) can be determined 
if we use the generalized random phase approximation 
for X (q,u) (Eq.|). 

As the density is lowered G(q) develops a peak cen- 
tered on q ~ 2Akp. When the height of the peak ex- 
ceeds unity there is the possibility that the denominator 
in the expression for x{q) (Eq. will pass through zero 
and in that case x(l) would diverge. Such a divergence 
would probably indicate a charge density wave (CDW) 
]15||. However we found that the peak in G{q) is not suf- 
ficiently high to offset the polarization which decreases 
rapidly with q, and we detected no CDW divergence in 

In Ref. M the local field factor was approximated with 
parameterized forms for parallel and antiparallel spin 
which had the correct large q asymptotic limit for G(q) 
fiof . We call this the Iwamoto Local Field Approxima- 
tion (ILFA). The parameters are expressed in terms of 
compressibility and spin-susceptibility sum rules in the 
long wave length limit. This approximation interpolates 
G(q) for intermediate values of q. For superconducting 
pairing the most important correlations are those acting 
over distances determined by the Fermi wave vector kp. 
Thus it is essential to determine the correlations in this 
intermediate q range accurately. We compare our results 
in the zero defect case with those obtained within the 
ILFA. 

We consider electron densities such that only the lowest 
energy sub-band is occupied. The effect of correlations 
on the effective electron-electron interaction 
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V(q)[l + V(q) X (q,<»)] (3) 



structure of the attractive region of V e ff(q,ep) is sensi- 
tive to correlations, as can be seen from the comparison 
with the V e ff(q,e F ) calculated within RPA (that is by 
setting G(q) = 0). V e ff(q,e F ) in RPA is only attrac- 
tive for q/k F < 0.2. This attraction is caused by the 
plasmon resonance and is not relevant here. When corre- 
lations are included, an attractive region develops over a 
wide range of wave-vectors. One already sees this effect 
within the ILFA (also shown), but the attraction there is 
much weaker than with our full G{q). The deep well near 
q/kp ~ 2.4 is the result of the peak in G(q) exceeding 
unity. This makes the term {V(q)[l - G(q)}} in Eq. | 
attractive. 




is shown in Fig. 1. We plot Re V e ff(q, oj=6f) for elec- 
tron density corresponding to r s = 10. The depth and 



FIG. 1. Effective electron-electron interaction V e ff(q, £f) 
(solid line) for zero disorder 7 = 0. Electron density corre- 
sponds to r s = 10. Dashed lines show V e //(q, ef) within RPA 
(no electron correlations) and the ILFA (see text), as labeled. 



In Fig. 2 we plot Re V e f f(q,Lu=e F ) for zero disor- 
der at various electron densities. Increasing the den- 
sity decreases the correlations, and we see that the at- 
tractive part of V e ff(q, e F ) weakens, especially around 
q/kp ~ 2.4. 

In Fig. 3 we show the effect of disorder. Re 
V e ff(q,uJ—ep) is plotted for electron densities corre- 
sponding to r s — 10 and 5 at two values of the disorder 
parameter 7 = and 7 = 2e F . This covers both the Si 
MOSFET samples where typically 7 = 2e F , and also the 
high mobility p-SiGc samples where the value of 7 <C e F . 
In the Figure we see that for r s — 10 the plasmon reso- 
nance around q/kp ~ 0.3 has been completely damped 
out by the time we reach 7 = 2e F . Nevertheless the at- 
tractive well for q/kp > 2 is not dramatically changed at 
this level of disorder. In contrast, at the higher electron 
density corresponding to r s — 5 the disorder weakens the 
attractive potential. 
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FIG. 2. Effective electron-electron interaction V e ff(q,ep) 
for zero disorder 7 = 0. Curves are for different r s as labeled. 



We have calculated the critical superconducting tem- 
perature T c using the Gor'kov equation for the gap func- 
tion. In the Cooper channel the irreducible interaction 
between two electrons at frequency u> and is given by 
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;V eff (q,in). (4) 



Here, the wave vector q is given in terms of the fre- 
quency u) as q = \/ p 2 + p' 2 — 2pp' cos 8, where p — 
y/2m*\uj + ep \ and p' = kp. The angle between p and p' 
is 8. Under the weak coupling approximation the critical 
temperature T c is given as 



T c = 1.134e F exp 



(1 + (^)) S 
(F 2 ) - A 



(5) 



The average (F) = doj(F(v) - F{0))/2u. We ne- 
glect any effect of disorder on G(q) so Xq=(F(0)) is inde- 
pendent of disorder. The numerical simulation data |l4| ] 
we have used is for disorder free systems. If one used a 
self-consistent calculation |ll| in which the mutual de- 
pendence of the disorder and the correlations were calcu- 
lated then the Ao would depend on disorder through the 
G(q). 

In Table 1 we show T c calculated for a range of r s and 
the levels of disorder 7 = and 2ep. Disorder decreases 
T c but the dependence is only significant at the higher 
densities: by r s > 10 T c is insensitive to this disorder. 
As r s decreases there is an increase in T c until r s = 6—7. 
Around this density T c passes through a weak maximum 
and then starts to decrease. 




FIG. 3. a. Effective electron-electron interaction 
V e ff(q, £f) for r 3 = 10. The curves are for two different val- 
ues of the disorder parameter 7/2eF as labeled, 
b. Same as (a) but for r s = 5. 





Tc (K) 


Tc (K) 




7 = 


7 = 2e F 


5 


0.62 


0.25 


6 


0.63 


0.43 


7 


0.58 


0.47 


8 


0.52 


0.46 


9 


0.46 


0.42 


10 


0.40 


0.38 



Table 1. Superconducting transition temperature T c as 
a function of r s for disorder levels 7 = and 2e F (see 
text). 

The discrepancy between our T c and Tl LFA calculated 
within the ILFA is much larger at r s = 5 than at r s = 10. 
At r s = 5 and 7 = ILFA gives T^ LFA = 0.27 K which 
is a factor of two smaller than our T c , while at r s = 10 
the Tl LFA = 0.34 K. The maximum in Tl LFA occurs 
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near r s — 8. 

In conclusion we find that electron correlations and the 
disorder typical of that found in the current experimental 
samples has a significant effect on the effective attractive 
interaction which acts between the correlated electrons at 
low densities. Interestingly while the quite high levels of 
disorder typical of those in the Si MOSFET experiments 
are strong enough to completely damp out the plasmon, 
the attractive interaction for q/kp > 0.5 persists because 
of the static correlations. The effective interaction is sen- 
sitive to the details of the correlations and our results in- 
cluding the full correlations significantly differ from the 
ILFA. The estimated superconducting transition temper- 
ature is sensitive to the details of the correlations and, 
at higher densities, to the disorder. 
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